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Abstract 

We discuss the issue of unitarity in particular quantum cosmological 
models with scalar field. The time variable is recovered, in this context, 
by using the Schutz’s formalism for a radiative fluid. Two cases are con¬ 
sidered: a phantom scalar field and an ordinary scalar field. For the first 
case, it is shown that the evolution is unitary provided a convenient factor 
ordering and inner product measure are chosen; the same happens for the 
ordinary scalar field, except for some special cases for which the Hamil¬ 
tonian is not self-adjoint but admits a self-adjoint extension. In all cases, 
even for those cases not exhibiting unitary evolution, the formal compu¬ 
tation of the expectation value of the scale factor indicates a non-singular 
bounce. The importance of the unitary evolution in quantum cosmology 
is briefly discussed. 
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1 Introduction 


Quantum cosmology faces difficulties that range from the technical to the con¬ 
ceptual point of view. It is usually based on the ADM formulation of the 
General Relativity theory m eh ®, leading to the Wlieeler-de Witt equa¬ 
tion, a functional equation defined in the superspace, the space of all possi¬ 
ble three-dimensional metrics on the space-like hypersurface which foliates the 
four-dimensional space-time. Technically, the foremost problem is to determine 
solutions of the Wheeler-de Witt equation in its complete form, since it is a 
functional equation with an infinite number of degrees of freedom. This sug¬ 
gests to freeze out an infinite number of such degrees of freedom, reducing a 
system to a few numbers of variables, leading to the mini-superspace setting. 
This drastic reduction implies that we ignore, for example, the violation of the 
uncertainty principle to those degrees of freedom that have been frozen out. On 
the other hand, the ADM formulation of General Relativity forms a constrained 
system, and the Hamiltonian constraints implies the absence of an explicit and 
clear time variable. On the conceptual side, moreover, we must face the problem 
of a unique system, the Universe at a whole, what leads to the inapplicability 
of the usual Copenhagen interpretation of quantum mechanics. 

The issue of the time variable has been addressed in many different manners. 
One can, for example, identify an internal time, determining the evolution of 
the system, through the deparametrisation procedure [5]. Or an external time 
(from the point of view of the gravitational system) may be introduced, for 
example, through the WKB approach (see, for example, reference [B] and refer¬ 
ences therein) or by considering matter fields. This last procedure has been used 
extensively in the literature, and one possibility is to consider a fluid with inter¬ 
nal degrees of freedom using Schutz’s description 0 0. This approach leads to a 
Schrodinger-like equation, since the momentum associated to the fluid variables 
appears linearly in the resulting Hamiltonian [5]. 

The interpretation problem is more delicate. This very important issue will 
not be treated in a direct way in the present text, but the reader can address 
himself to many text about this subject in the literature. See, as example, 
references 0 uni nn nil [inj and references therein. 

Our goal in this paper is to investigate another difficulty that appears even 
when the mini-superspace approach is employed and the time variable is iden¬ 
tified through, for example, the Schutz’s variable: the unitary evolution of re¬ 
sulting quantum system. It has been shown that when just the Einstein-Hilbert 
Lagrangian and a perfect fluid matter component is considered, the resulting 
Hamiltonian can be made self-adjoint (unitarity is assured) under the hypothe¬ 
sis of a maximally symmetric spatial section M- However, when anisotropy is 
taken into account, in the same context, the Hamiltonian is generally no more 
self-adjoint m■ However, recently, it has been shown that a convenient choice 
of the ordering factor associated with the gravitational operators may restore 
the self-adjoint character of the Hamiltonian [16]. 

The breakdown of the unitary evolution when the time variable is recovered 
through the Schutz variable occurs also when scalar fields are considered to- 
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gether with the Einstein-Hilbert Lagrangian, even in the isotropic e homogenous 
case mm- Curiously, in this case, and using the Schutz formalism to recover a 
Schrodinger-like equation, a phantom scalar field (with negative kinetic energy) 
assures the positivity of the total energy, while for a normal scalar field (with 
positive kinetic energy) this positivity of the energy is not assured, since the 
Hamiltonian has a hyperbolic signature for this case. It will be shown, however, 
that a convenient factor ordering (similarly to what has been discussed in the 
anisotropic case, reference m) together with a convenient measure associated 
with the inner product guarantees a self-adjoint Hamiltonian (hence a unitary 
evolution) for a phantom scalar field. The same occurs for the normal scalar 
field, but for special cases (connected with the ordering factor) the Hamiltonian 
is not self-adjoint but admits a self-adjoint extension. We will not consider here 
the issue of unitary evolution for other approaches to recover the time variable 
as, for example, by using the scalar field itself as the time coordinate m- 

This paper is organised as follows. In next section with consider a scalar- 
tensor theory - specifically, the Brans-Dicke theory - coupled with a radiative 
fluid. By performing conformal transformation, the theory is re-written in the 
so-called Einstein’s frame. In section 3, solutions for the Schrodinger-like equa¬ 
tion are obtained, and it is shown that they do not display a unitary evolution. 
Even though, formal predictions for the evolution of the Universe are obtained 
revealing a non-singular behaviour. In section 5 it is shown how to recover 
unitarity at least for the anomalous, phantom configuration of the scalar field. 
The general self-adjoint issue is considered in section 6. In section 7 we dis¬ 
cuss the self-adjoint extension for the particular cases where the Hamiltonian 
is not self-adjoint but admits such an extension. In section 8 we present our 
conclusions. 


2 Scalar-tensor model with radiative fluid 

The prototype of the scalar-tensor gravity formulation is the Brans-Dicke theory 
[Mj, represented by the Lagrangian, 

£ = — |+-£m, (1) 

where L m is the matter Lagrangian. We will suppose that the matter Lagrangian 
refers to a radiative fluid, having conformal symmetry. This Lagrangian includes 
as a particular case the string dilatonic Lagrangian for which uj = — 1 [2Tj . 

The Lagrangian m defines the theory in the Jordan’s frame. Performing a 
conformal transformation such that g = <(> -1 ( 7 ^„, we transpose the action UJ 
to the corresponding expression written in the Einstein’s frame |22j : 

U = ^| J R- W ^J+£ m , (2) 

where u> = Co + 3/2. Remark that the matter Lagrangian is not affected by the 
conformal transformation since we are considering a radiative fluid. 
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We will consider from now on that the line element in the Einstein’s frame 
can be written as, 

ds 2 = N(t) 2 dt 2 — a(t) 2 'yijdx l dx - 7 (3) 

where N(t) is the lapse function, a(t) is the scale factor and 7 y is the induced 
metric of the homogeneous and isotropic spatial hypersurfaces with curvature 
k = 0, ±1. For simplicity, we will fix k = 0. With this metric, the gravitational 
Lagrangian becomes, 


£g 


boa 3 

N 


-6 



aN 1 </> 2 'I 

a N (f> 2 j 


(4) 


where Vq is a constant and can be interpreted as the physical volume of the 
compact universe (in this case a three-torus) divided by a 3 . Since there is an 
identical multiplicative constant in front of the matter Lagrangian we can drop 
it from our analysis. Discarding a surface term, the gravitational Lagrangian 
can be written as, 

£ G = ^{ 6 aa 2 -w« 3 ^}. (5) 

Defining, 

cr =V / MW, ( 6 ) 

we obtain, 

-Cg =-^-{saa 2 - ea 3 d 2 }, (7) 


where e = ±1 according u is positive (upper sign) or negative (lower sign). 
The canonical momenta associated with the scale factor and the scalar field are 
respectively: 


aa 

Pa = 12 — 

Pa N 


= - 2t ~ 


leading to the following expression in terms of the conjugate momentum: 


£G = Paa + PaV - N 


1 Pi Pl\ 
24 a 4a 3 J 


( 8 ) 


(9) 


Considering a radiative matter component (for the computation of the con¬ 
jugate momentum associated with the fluid, see references 0 fl4l). the total 
Hamiltonian is: 


24 a 4 a 3 a 


The resulting Sclirodinger-like equation is, 

d 2 ^ e d 2 ^ 


da 2 

where we made the redefinition 


da 2 


= i 


.d^ 
dT ’ 


a and ^ —> T. 


( 10 ) 


( 11 ) 
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3 Cosmological scenarios 


In the reference [18], the quantum cosmological model defined by the Schrodinger- 
type equation m has been studied. We will review in this section the procedure 
and results obtained in reference [IS] . 

In order to obtain treatable expressions, the ordering ambiguity of operators 
has been exploited, and equation m has been rewritten as, 


d 2 T 1 , 

—— d„A’ 

oa- a 




e d 2 ^ 
a 2 da 2 dT' 


( 12 ) 


For e = — 1, a condition that assures the positivity of energy, equation m 
admits a solution in terms of stationary states of energy E: 


v]/ = AJ v (\fEa)e 


ikcr —iET 

& i 


= k 


(13) 


with A being a normalization constant and k is a separation constant. 

A particular wavepacket can be obtained by a convenient superposition of 
the constants E and k , as it is described in [18]. The final result is 


T(a,CT,T) = C 


a 2 

e ~4 B(T) 

B{T) g a (a, B.a) 1 


(14) 


where 


B(T) = ( 1 + iT), 


g a {a,a,T) =-a + In 


2 B(T) 


± 1 ( 7 , 


(15) 


7 and a being constants connected with the gaussian-type superposition, and 
C is a normalisation constant. 

The norm of the wavefunction m can be calculated explicitly: 


N 


C 2 


10 


oo /»+oo 

-oc *** dada= {BB*y/ 2 ’' J 0 


o~lU 


rTT 


a + In ( f 


du 


C 2 


(B B*) 1 / 2 


7T/l, 


(16) 


where h is the definite integral, 


h 



(17) 


The norm is clearly time-dependent: the corresponding quantum model is not 
unitary. Even though, the expectation value for the scalar field can be formally 
computed, leading to the expression, 

< a >oc (q 2 ±T 2 ) 1/2 . (18) 
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The same result is, essentially, obtained through the computation of the bohmian 
trajectories. The expectation value for the scale factor indicates a non-singular 
bounce. As it is shown in reference [IS], the expectation value of the scalar field 
a is time-dependent, reading 

< cr >oc arctan ^ . (19) 


4 Recovering unitarity 


Let us suppress the ordering factor introduced previously. The Schrodinger 
equation is 


- + ^2^1' = id T ’F, 

with e = ±1. For a stationary state, 


if = cj)e- iET , 


leading to 

- dl<P + ^0 = E<p. 

Using the separation variables method, such that 
0(a,cr) = X(a)Y(a), 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


the equation becomes, 


, e Y 

+ ~ E ' 


(24) 


where the primes mean derivative with respect to a and the dots mean derivative 
with respect to a. This equation can be rewritten as, 


ey = a'\E 


A" 

~Y 


= — ek 2 , 


(25) 


where k is a constant of separation. In this case, the function Y satisfies the 
equation, 

Y + k 2 Y = 0, (26) 

with the solution, 

Y = Y 0 e ika . (27) 

The equation for X reads, 

X"+^E + e^PjX = 0. (28) 
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The soluton is, 


X(a) = y/aJ u (V~Ea), v = ^ — ek 2 . 


(29) 


The total wavefunction is, 

= m^J u (y/Ea)e-^ ka+ET) . (30) 

Let us construct the wavepacket. First we write x = \[~E. Then, we can 
made the following superposition: 

/ +oo r + oo 

/ A{k)x v+1 e~ Bx J v {xa)e~ lktT dkdx, (31) 

- oo J 0 


where, as before, 


B = 7 + iT, 


(32) 


7 being a positive constant (in order to assure the convergence of the integral), 
playing the same role as the previous section. The function A(k) will remain 
for the moment unspecified, but it must decay exponentially for large k in order 
also to assure the convergence of the wavepacket. The integral in x can be 
performed leading to '231 . 


f + OO 




A(k)e 


—ika 


a ' 


H-l/2 


( 22 ?) 


v+i 


exp ( — 4g I dk- 


(33) 


From now on, let us consider the case e = — 1, such that v is real. The 
modulus of the wavefunction reads, 


/ +oo r+oo t 

/ A(k)A*(k')e-^ k - k '' ),T B- , '- 1 B*~ v ~ 1 a v+v ' +1 

-oo J — oo 




x exp 


7GT 
22? 


dkdk ', (34) 


where B = BB* = y 2 + T 2 . 

The integration in a in the interval —oo < a < +oo implies a delta function, 
S(k — k'). After integrating in k', we find 


/*+oo 


p + OO 




\A(k)\ 2 B- 2v - 1 a 2v+1 e~*B dk. 


(35) 


Performing the change of variable, 


V = 


Vb' 


we can write the norm as, 

C + OO /* + oo 


r too /‘too 

N= / \A(k)\ 2 e~^ v y 2v+1 dydk. 

J 0 J —oo 


(36) 
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The norm is time-independent. After integrating in y. it can also be written as, 


N = \ y_ + °° \A{k)\ 2 r(^-^j dk. (37) 

Due to the asymptotic properties of the T’s function, the integral converges if 
A(k) behaves, for example, as e~ k for large values of k. 

The expectation value for the scalar field is given by the expression, 

^ p- 1-00 p+oo 

<a>=— / / da da. (38) 

N J- oo Jo 

Using the same procedure as before, we end up with, 



which is essentially the same prediction as in the previous non-unitary case. 

We can also evaluate the expectation value for the scalar field. In fact, to 
obtain the expectation value of a we must compute 

nOO /» + 00 

< tr > oc / / ^(a, cr)*<7'h(a, a)da da. (40) 

J 0 J— oo 

Using (1551) . we can write 

r+°° a v+1/ 2 / a 2 \ 

a ' i ’ = i J-oo A ^ dke 1 J ] (2 BY+ 1 exp (~ I BJ dk ' (41) 


which can be integrated by parts. Using this expression, the integration in a 
leads again to delta function. Redefining, as before, the integration in a, we 
obtain the following expression, 


< (T >T 



A k {k) 

m 

2^+1 

dy dk , 


(42) 


where the subscript k indicates derivative with respect to this parameter. After 
expressing B in polar form, we obtain, 


< a >t= + o\ arctan 



(43) 


where erg and <j\ are constants. This is essentially the same result obtained 
in the previous section for the expectation value of a. Remark, however, that 
< a >= 0 if A(k) is an even function, such that A(k) = A(—k). 
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5 Self-adjoint ness 

The results shown before reveal that the self-adjointness property of the effec¬ 
tive Hamiltonian may depend on the ordering factor. The overall situation is, 
however, more involved. The Hamiltonian used in equation (1121) is, 

H = -dl - ±0 O + (44) 

When e = — 1, this operator looks like that one of the two dimensional problem 
written in polar coordinates. In fact, if we make the identification a = r, a = 6, 
considering the variable 9 as periodic, such that we can define the cartesian 
coordinates x = rcosO and y = rsind, the Hamiltonian takes the form, 

H = -dl- d 2 y , (45) 

with —oo <x,y < +oo, which is clearly self-adjoint. 

What makes the Hamiltonian (HU) not self-adjoint? We could think that is 
the fact that the coordinate a there is not periodic. But, we remark that in 
passing from cartesian to polar coordinates, there is a Jacobian factor that it 
has not been used in computing the inner product in the Hilbert space for the 
Hamiltonian (|44l) . The clue to this problem seems to be in the ordering factor 
and in the measure of the inner product, as point out in reference m in the 
context of anisotropic models. 

Let us be more specific. Let us consider the Hamiltonian with a ordering 
factor labeled by q: 

H = -dl - *d a + \dl, (46) 

a or 

where we have restored the factor e = ±1. This Hamiltonian is symmetric under 
the inner product defined by, 

poo p-\-oo 

(0,4')= / / 0*4' a q dada, (47) 

J 0 J —oo 

if the functions 0 and 4', as well as their first derivatives, are null in the extreme 
of the interval. 

To be symmetric is a necessary but not sufficient condition for the operator 
H to be self-adjoint: the domain of the operator and of is hermitian conjugate 
must be also the same. The self-adjoint character of a symmetric Hamiltonian 
can be obtained through the deficiencies indices of von Neumann [241 (251 . 

The deficiencies indices are defined as follows. Consider the eigenvalue prob¬ 
lem, 

H 4- = ±i4>. (48) 

Let us call n+ and n_ the number of square integrable solutions of the eigenvalue 
problem given by equation (1481) when the upper and lower sign in the right 


hand side is chosen, respectively. If n + = n_ = 0, the operator H is already 
self-adjoint; if n + = n- ^ 0, the operator is not self-adjoint but it admits self- 
adjoint extensions given by some restrictions in the wavefunctions; if n+ ^ Ti¬ 
the operator is not self-adjoint and it does not admit any self-ajoint extension. 
Hence, we must solve the eigenvalue equation, 

- dl'i> - -da'b + Ar'P = ±H1t = 77 ^, (49) 

a a z 

where we have defined r] = ±i. The solutions of this equation can be written as 


= a p 


HPiy/rja) + H?\Vva) 


_=h ika 


(50) 


where, 

P= v = Vp 2 - ( 5 1 ) 

k is a separation constant which must be real (otherwise, in the corresponding 
Schrodinger equation, we would have just non-integrable solutions, with diver¬ 
gent norm), and H^' 2 \x) are Hankel’s functions of first and second kind. 

Let us write the functions, 

= a p (r)a)e ika , (52) 

¥ 2) = a p H^(r/a)e ik<T . (53) 


These solutions contain the separation parameter fc, which parametrises the 
plane wave behaviour in terms of the coordinate er. In order to give a physical 
meaning to such structure is necessary to construct a wave packet related to 
these plane wave solutions. This procedure has the purpose to avoid unphysical 
behaviour, but it has no major impact on the considerations on the self-adjoint 
character of the operator. Hence, the solutions take the form, 


/ +oo 

A(k) a p Hl 1 \rja)e ik ' 7 dk , 

-OO 

/ +oo 

A(k) aP (r]a)e ikcr dk , 

-OO 


(54) 

(55) 


where A{k ), as before, is a superposition factor which satisfies the condition to 
go to zero sufficiently fast at both infinities. 

We must investigate the norm of the wave functions (1541551) . given by 

/ +oo /»+oo + 

/ (a, cr)(a, a) a q da da. (56) 

-oo J 0 

First, let us consider the integration in a. We have: 

/*+oo 

I* = *£ 2) '*£ 2) da 


9 








p+oo /*+oo /»+oo 


J— oo J — oo J—oo 

/ +oo />+oo 

-oo —oo 
r+oo 


A(k)A* ( k)a 2p i# 2) ( na) h£’ 2) (pa) e i(k ~ k)a dk dk da 
A(k)A* (k)a 2p ( pa(pa ) 8(k — k) dk dk 

/ A~00 

\A(k)\ 2 a 2p (rja)Hj) 1,2 \rja) dk. (57) 

-oo 

The norm becomes, 

OO /»-l-oc 

| A(k)\ 2 a 2p+q H^’ 2) *(pa)Hj > 1 ’ 2 \pa) dk da 

2n / \A(k)\ 2 a hA’ 2 ' 1 (pa)H^ 1,2 \pa) dk da. (58) 

J 0 J —oo 

We may investigate separately the norm in two limits, when a —> 0 and 
a —> oo. 


/ + CXD r+OO 

/ 

-oo J 0 

/*+oo /»+oo 

= 2, / 

</ 0 J —oo 
/»+oo /»+oo 


• a — > 0. In this case, 

H^irja) ~ a~ v . ( 59 ) 

The norm behaves as, 

/ +oo 

\A{k)\ 2 a 2[1 -AWdk., (60) 

-OO 

where we have clearly indicated that v depends on k. We have two different 
situations. 

1. If e = —1, v € 3? + . In fact, \p\ < v < oo. There is no square- 
integrable eigenfunction in this case. This indicates that the Hamil¬ 
tonian operator is self-adjoint. 

2. If e = 1, v = \Jp 2 — k 2 , and the divergence of the norm in this limit 
depends on the ordering factor and on the interval of k. Suppos¬ 
ing, for example, that the superposition factor A(k) is defined in all 
interval of k, there is always divergences unless —y/l — p 2 < k < 
+vl —p 2 . In order to assure the existence of this interval of values, 
it is required that p < 1. 

• a —> oo. In this limit, 




TyAja 




(61) 
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where the upper (lower sign) in the exponential corresponds to the Han- 
kel’s function of the first (second) kind. In this limit, the norm of the wave 
functions takes the form, 

9 r+oo r+oo 

iV±~-/ / \A(k)\ 2 ae^^^-^-^dkda, 

77 Jo J — oo 

where the symbol (±) corresponds to the two types of eigenfunctions of 
the equation (05]). It is clear that for each type of Hankel’s function there 
is one eigenfunction that converges (leading to a finite contribution to the 
norm) and the other that diverges (leading to a infinite contribution to 
the norm). 

Hence, with these informations for the two asymptotic regimes, we can con¬ 
clude the following. 

1. For e = — 1 we have just divergent eigenfunctions corresponding to the 
equation (1481) due to the behaviour in the limit a —> 0. This implies the 
deficiency indices are such that n. (_ = n_ = 0. The Hamiltonian operator 
is self-adjoint. 

2. For e = 1, there are two possibilities. 

(a) For an ordering factor such that p > 1, we find, as previously, diver¬ 
gent eigenfunctions for equation (1481) implying n_|_ = n_ = 0 and the 
Hamiltonian operator is already self-adjoint. 

(b) For an ordering factor such that p < 1, the wave functions are con¬ 
vergent in the limit a —> 0, but in the limit a — > oo we have one 
divergent and one convergent function for each sign in the equation 
(1481) . These asymptotic behaviours imply n. (_ = n_ = 1. Hence, the 
Hamiltonian is not self-adjoint but admits a self-adjoint extension. 

6 Self-adjoint extension 

The conclusion of the preceding section was that the case e = — 1 leads to a 
Hamiltonian that is already self-adjoint, but for the case e = 1 the Hamiltonian 
is self-adjoint if the factoring order p is greater than 1, while for p < 1 the 
Hamiltonian is not self-adjoint, admitting however a self-adjoint extension. The 
goal of the present section is to verify the conditions on the wave functions in 
order to have this self-adjoint extension when e = 1 and p < 1. 

The self-adjoint extention are given by the functions (p such that, 

wm = (Hm, ( 62 ) 

where = ’F+ + for some A. We will keep, just for generality, the value of 
e undefined. We have then, 

/ oo poo 

/ ip* (H4>) a q dada 

-oo J 0 


li 


= J J ip* [-d%<p- ~d a <l>+ a q dada. (63) 

Let’s compute the integral each term separately. The first terms read, 

j J -ip* (d%<p) a q dada = J [-ip* a q d a <p\ a a ..Z™ - [~d a (ip*a q ) </>]“ da 

—d 2 ( ip*a q ) (pdada 

I {[(dalp*) (p-ip* {da<P)\ a q + qa q ~ 1 ip*(p} a a ““ da 

J — OO 

+ f f — (d 2 ip*) (pa 9 — 2q (d a ip*) (pa 9-1 — q (q — 1) ip* (pa q ~ 2 dcida. (64) 


pOO p OO 


J — oo J 0 

p OO 


The second term is given by, 



a q dada 



[-qr<Pa q -Xl7 da 


n oo 

q (d a ip*) (pa q ~ x 


+ q (q — 1) ip* (pa 9 2 dada, 


(65) 


Finally, the third term is, 

/ [ ip* (-^9 2 cp) a q dada = f [ip* (d a <p) - (d a ip*) 4>ZZ-oo a<?da 

J-ooJO va ' JO a 

/ oo /»oo 

/ “ 2 <t>a q dad(j. (66) 

-oo Vo ^ 


Putting it all those expressions together, we obtain 

/ oo 

{[(da-ip*) cp-ip* ( da(p)\ da 

-OO 

+ / 00 4[^(M-(a^*)^! 0 00 a , do 

Jo a 

+ J j (-d^ip* - ^ d a ip* + Ajdpip*^ (pa q dada. (67) 


Then, in order to get (ip\H<p) = (Hip\cp), the function (p must satisfy the bound¬ 
ary conditions: 


{[(d a ip*) (p - ip* {d a (p)} a q } a a JZ = 0; (68) 

[V>*(^)-(d^*)C=-oo = 0. (69) 

Nevertheless, since (p is already square integrable, and goes to zero sufficiently 
quickly in the extremes of the interval of a , and the same happens for a —> oo, 
the only condition that is not automatically satisfied is 

lirn {[(d a ip*) <p-ip* (d a <P)\ a 9 } = 0. (70) 

a— 
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Remembering that we found previously, 


*+ = 

= 


A (k) aPH^ (Via) e ifc<T d/e 

/ OO 

A{k)aPHl 2) (V^ia) e ika dk , 

-OO ^ 7 


-OO 

/•OO 


then, 


= 

V>* = 


' —OO 
/•OO 


A (A) a p (Via) + A 
A* (jfe) aP HW (Via) + Aif< 2) 


S ifeff d/c; (71) 


e~ ika dk. (72) 


The condition is given for a —>• 0 and, in this limite we have the expression 
(fSiip . Therefore 

/ OO 

A* (jfe) aP (a~ v + A a~ v ) e~ lka dk. 

-OO 

Hence, we must have, 

‘ /*oo 

/ A* (jfe) aP (a~ v + Aa - ") e~ lka dk (j)a q 

—OO 
■ /* OO 

/ H* (jfe) a p (a - " + Aa - ") e~ ikG dk (, 9 a </>) a 9 

_«/ —OO 


lira <9 a 

ct—^0 


— lim 

ci—>-0 


= lim 

a—>-0 


— lim 

a—>-0 


/ oo 

H* (jfe) (p - i/) aP-'-ie-^ + H* (k) A (p - i/) aP- v ~ 1 e- ikG dk 

-oo 

/ OO 

H* (jfe) a p (a - " + Aa"") e -ife<7 djfc (d a </>) a q 

-oo 


^>a 9 


= lim 

a—>-0 


ip-v) 


/ OO 

H* (jfe) a p (a"" + Aa" 17 ) e~ ika dk 

-OO 


(j) a 9 


— lim 

a—>-0 


= lim 

a—>-0 


/ oo 

H* (jfe) a p (a - " + Aa -1 ') e -lfeCT dfc 

-oo 

/ OO 

H* (jfe) a p+9 (a -1 ' + Aa-") e~ lkG dk 

-OO 


(5 a </>) a 9 

ip ~ V <t> 


- ( d a (f >) 


(73) 


Since p + q = p + (—2p — 1) = — p — 1 < 0 and v > 0, the integral term diverges 
and the condition reduces to 


lim 

Ob - >0 


(p-v)<t> 


- {da4>) 


= 0. 


(74) 


Then, the self-adjoint extension is given in the subset of the square integrable 
function that satisfies the condition (1711) . 
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7 Conclusions 


In this paper we have addressed the problem of unitary evolution of isotropic 
and homogenous quantum cosmological model in scalar-tensor theories of grav¬ 
ity, with a time variable defined by the matter fluid. As it happens in the case 
of the anisotropic universe without scalar field m, the unitarity is assured if an 
ordering factor and an inner product measure are conveniently chosen. We have 
considered the case where the matter fluid is radiation, hence, bearing conformal 
invariance. This is somehow crucial in our analysis since we have consider ini¬ 
tially a scalar-tensor theory defined in the Jordan’s frame, and transformed it to 
a minimal coupling configuration, the Einstein’s frame. Hence, our conclusions 
are restricted to such framework. 

The analysis we have performed uses the traditional techniques of quantum 
mechanics, including the computation of the deficiency indices to determine the 
self-adjoint character of the Hamiltonian, as well as the verification of the self- 
adjoint extension for the case the Hamiltonian operator is not self-adjoint but 
admits a self-adjoint extension [24 ] [25 ] . 

There is an essential point we have touched only superficially in the present 
work: given a self-adjoint operator, that assures a unitary evolution for the 
system, what are the predictions to evolution of the Universe? In reference [18] , 
where a non-unitary model has been considered, the predictions for the evolution 
of the universe were obtained by considering the bohmian’s trajectories, and 
a non-singular universe, with a bounce, has been obtained. After restoring 
unitarity, at least for a particular case, essentially the same predictions were 
obtained either for the scale factor and for the scalar field. 

The results described above may lead to a different issues: what is the role 
of unitary evolution in the universe and how it is connected with the choice of 
given time coordinate. In references d 1271, for example, scalar-tensor theories 
were analysed without introducing an external time variable connected with a 
dynamical fluid, and the predictions obtained through the WKB method and by 
bohmian trajectories indicated a non-trivial behaviour for the scalar field and 
for the scale factor. It is important to compare these different procedures. 

The role played by the unitary evolution, which supposes the choice of a 
variable parametrising the evolution of the system, and the equivalence be¬ 
tween different choices for this time variable, must be discussed more deeply. 
The consequences of not having a positive defined energy when the scalar field 
has positive kinetic energy (but with a self-adjoint Hamiltonian) must be also 
investigated. We hope to address these questions in future works. 
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